Quasiparticle excitation in and around the vortex core of underdoped YBa 2 Cu 4 08 
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We report a site-selective 17 O spin-lattice relaxation rate Tj~ in the vortex state of underdoped 
YBa2Cu4 08. We found that T^ 1 at the planar sites exhibits an unusual nonmonotonic NMR 
frequency dependence. In the region well outside the vortex core, T^ 1 cannot be simply explained 
by the density of states of the Doppler-shifted quasiparticles in the d-wave superconductor. Based 
on T-^ 1 in the vortex core region, we establish strong evidence that the local density of states within 
the vortex core is strongly reduced. 



The electronic structure in the vortex state is a funda- 
mental problem in the physics of superconductors with 
unconventional symmetries. Recent experimental results 
have revealed that the quasiparticle (QP) excitations in 
the vortex state of high-T c cuprates are very different 
from those of conventional s-wave superconductors. De- 
spite extensive studies on this subject, the question of the 
QP excitations both inside and outside the vortex core 
is still far from being settled. In fact, STM experiments 
have revealed many unexpected properties in the spec- 
trum of the vortex core jjj. For example, a zero energy 
peak in the tunnel conductance due to the extended zero 
energy QP states, which is expected in the conventional 
theory of the d-wave superconductor, is not observed at 
the center of the vortex core. Although the novel elec- 
tronic structure of the vortex core has been discussed in 
terms of several intriguing models, it remains controver- 
sial [§-f§. O n the other hand, well outside the vortex 
core, the QPs are considered to be in extended states 
due to the presence of gap nodes. The dominant effect 
of these delocalized states in the magnetic field is the 
Doppler shift of the QP energy spectrum, which occurs 
due to the presence of a supercurrent flow around each 
vortex and generates a nonzero density of states (DOS) 
at the Fermi surface ||. The H- and T-dependencies of 
the specific heat 0, thermal conductivity and NMR 
Knight shift I| at low temperatures have generally been 
accepted as strong evidence for this effect. However, 
rather surprisingly, neither the local DOS with a 1/r 
asymptotic tail expected from the Doppler shift nor the 
4-fold symmetry of the extended QPs has been observed 
in the STM spectrum GJ]. Thus, a discrepancy still exists 
between the STM findings and the results of other ex- 
periments including thermodynamic measurement. The 
situation therefore strongly confronts us with the need 
for a powerful probe of the QP excitation with spatial 
resolution. 

The measurements of the NMR spin-lattice relaxation 
rate TV" 1 provide microscopic information on the low en- 
ergy QP excitations in the superconducting state. Re- 



cent experimental and theoretical studies have estab- 
lished that the frequency dependence of T^ 1 measure- 
ments serves as a site selective probe for the electronic 
spectrum in the vortex state of the high-T c cuprates 
||lo| fl3|| . Although the site-selective NMR, compared to 
STM, does not possess an atomic resolution , it has two 
advantages. First, it is a bulk probe, free from surface ef- 
fects. Second, as pointed out by several authors H,E5[, 
the tunneling conductance obtained by STM is not sim- 
ply related to the local DOS owing to the very anisotropic 
tunneling matrix element along the c-axis. The site se- 
lective NMR, on the other hand, especially at the oxygen 
sites in the plane where the antiferromagnetic spin fluctu- 
ation is cancelled, provides direct information on the local 
DOS |r|[L7|]. Very recently, the frequency dependence of 
T^ 1 at the 17 sites (denoted as 17 T 1 _1 ) in the vortex 
state of optimally doped YBa2Cu30y has been reported 
JTo| . A key feature of the data is that 17 T 1 _1 at the planar 
sites increases monotonically with NMR frequency well 
outside the vortex core, which is consistent with the de- 
localized QP DOS induced by the Doppler shift. In this 
Letter, we attempted to measure the site-selective NMR 
in underdoped YBa2Cu40s extending the measurement 
to the vortex core region which has not yet been reported. 
We have found an unusual NMR frequency dependence 
of 17 T 1 ~ 1 at the planar sites in the vortex lattice state, 
which is markedly different from that of YBa2Cu307. On 
the basis of these results, we discuss the QP structure in 
and around the vortex core of YBa2Cu40s. 

High quality underdoped YBa2Cu40g with T c =81 K 
has been synthesized under a high-oxygen-pressure at- 
mosphere. Previous NMR studies have shown that the 
normal state spin fluctuations at the planar Cu and O 
sites for YBa2Cu40s are similar in behavior to those 
for the 60 K phase of YBa 2 Cu 3 6+ z §J]. The NMR 
data were obtained on an aligned powder of underdoped 
YBa2Cu408, which was isotopically enriched with 17 0. 
The 17 spin echo signals were obtained by a conven- 
tional pulse spectrometer under a field cooling condition 
(FCC) in a constant field (Hq=9A T) applied parallel 
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to the c-axis by using a highly homogeneous supercon- 
ducting magnet. We stress here that the NMR measure- 
ments under the FCC are important because the Bean 
critical current associated with the sweeping magnetic 
field not only produces the field gradient in the crystal 
but also seriously influences 17 T 1 _1 by producing an ad- 
ditional Doppler shift to the QPs. 

The final NMR spectra were obtained by convolution 
of the respective FT-spectra of the spin echo signals mea- 
sured by a step of 50 kHz. Given that the spectra of the 
central transition of the chain 0(1), the planar 0(2,3), 
and the apical 0(4) overlap at low temperature, 17 T\ 
is determined by the recovery curves of the quadrupole 
split 0(2,3) satellite. The inset of Fig. 1(a) shows the pla- 
nar 17 (0(2) and 0(3)) spectrum of the satellite tran- 
sition at 100 K. The slightly-split between the satellite 
spectra of the oxygen sites of 0(2) and 0(3) is due to 
a small difference in electric field gradient. The proce- 
dure to determine T\ is as follows. The spin echo in- 
tensities are measured as a function of time after sat- 
uration pulses. Then the FT spectrum of each echo is 
obtained. Each frequency component of the FT spectra 
is fitted individually by R(t) = (M (oo)-M(i))/M (oo) = 
l/35exp(-2W*)+3/56exp(-6Wi)+l/40exp(-12Wt) + 
25/56 exp(-20Wi) + 25/56 exp(-30Wt) 

for the (3/2<->l/2) transition. Here, we calculated 17 Ti 
(=1/2W) for each frequency point at the 3 kHz interval 
with a gaussian weight function of er=3 kHz. 

The local field profile in the vortex state is obtained 
by approximating Hi oc (r) with the London result, 



Hi oc (r) = H ^2 CX P (~ iG ' r ) 



C M-e ab G 2 m 



(1) 



where G is a reciprocal vector of the vortex lattice, 
r | is the distance from the center of the core, ^ a b is 
the in-plane coherence length, and X a b is the in-plane 
penetration length. The thin solid lines in Figs. 1 (a) 
and (b) depict the histogram at particular local fields 
which are given by the local field distribution f(Hi oc ) = 
Jo S[Hi oc (r) — Hi oc ]d 2 r , where f2 is the magnetic unit cell. 
In the calculation we used £=15A and X a b = 1600 A re- 
ported in Ref. and assumed the square vortex lat- 
tice. The magnetic field is the lowest at the center of the 
vortex square lattice (C-point in the inset of Fig. 1(b)), 
and it is the highest at the center of the vortex core (A- 
point). The intensity of the histogram shows a peak at 
the field corresponding to the saddle point of the vortices 
(B-point). Because the local NMR frequency is propor- 
tional to the local magnetic field, we are able to produce 
a spatial image of the low-lying QP excitation spectrum. 

The solid lines in Figs. 1 (a) and (b) represent the 
NMR spectrum at the planar and the apical sites in the 
superconducting state, respectively. A clear asymmetric 
pattern, which originates from the local field distribution 
associated with the Abrikosov vortex lattice (the 
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FIG. 1. The planar 17 0(2, 3) (a) and the apical 17 0(4) 
(b) satellite (3/2<->l/2 transition) spectrum (solid lines) at 
T=20 K and 10 K, respectively. The open circles show the 
T{~ at corresponding sites. The thin solid lines represent 
Redfield patterns for planar 17 0(2), 17 0(3) and apical 17 0(4) 
sites. The broken lines represent the simulation spectrum 
convoluted with the Gaussian broadening. For details, see 
the text. The inset of Fig.l (a) shows the satellite spectrum 
of 17 0(2,3) at 100K. The inset of Fig.l (b) shows the image of 
field distribution in the square lattice of vortices. A-, B-, and 
C-points represent the vortex core center, the saddle point 
of the field distribution, and the middle of the vortex square 
lattice, respectively. 

Redfield's pattern), is observed well below T c in both 
sites. The broken lines in Figs. 1 (a) and (b) depict the 
spectrum convoluted with the Gaussian broadening func- 
tion, f(H loc ) = e^ 2H ioJ^ . (The broken line in Fig. 1(a) 
is obtained by the sum of the spectrum at the 0(2) and 
(3) sites.) Using cr=24 kHz, the theoretical curves well 
reproduce the data in both planar and apical sites. The 
open circles in Figs. 1 (a) and (b) display the NMR fre- 
quency dependence of 17 T± at the planar (denoted as 
17 T 1 ~p i ) and at the apical ( 17 T 1 ~ C J P ) sites, respectively. A 
clear nonmonotonic frequency dependence is observed at 
the planar sites. With increasing frequency, 17 T^~^ first 
decreases, shows a minimum near the B-point, then in- 
creases. Finally, at high frequency when approaching 
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FIG. 2. 1/Ti of 17 0(2, 3) at planar sites, 
tion of frequency for various temperatures 
the half width at half intensity of the spectrum of 1 '0(2) as 
a function of temperature. T c is the superconducting tran- 
sition temperature, and T me it is the vortex lattice melting 
temperature. [^o| 

the vortex core regime (A-point), 17 T' 1 _ p i again decreases. 
This frequency dependence is in contrast to the result 
of Ref . JhJ , in which 17 T , 1 ~^ shows monotonic increases 

with frequency. In contrast to 17 T^ V 17 ^i~ap snows little 
frequency dependence. 

Figure 2 shows the frequency dependence of 17 T^ pl at 
various temperatures. The inset shows the half width at 
half intensity of the spectrum of 17 0(2) as a function of 
temperature. One can clearly see a slight increase of the 
line width below T c and the large broadening due to the 
vortex lattice below T me i t ~ 45 K. (20) 17 t> ecomes 
nearly independent of frequency at ~ 45 K, above which 
the vortex lattice melts into disordered vortex liquid. The 
frequency dependences of 17 T^~ pl are remarkable below 



FIG. 3. T-dependences of 1/Ti of planar 17 0(2,3) and api- 
cal 17 0(4) sites, v (A) represent the largest (the smallest) 
values of of 17 0(2,3) across frequencies. The dotted 

line for the ir O(4) data represents the TiT=const. relation. 

smallest values is very similar to results reported pre- 
viously. p8[ Because T\ is determined from averaged re- 
covery curves over all frequency components by a conven- 
tional method, T\ is attributed to the value from the most 
dominant components of the signals, namely, in the su- 
perconducting region far from vortex cores (around the B 
point). With decreasing T, 17 decreases faster than 
T above T c due to the pseudogap, and does not show 
an anomaly at T c Jl(|. On the other hand, l/ 17 Ti, p is 
nearly proportional to T except at ~45 K, where a peak 
in 17 7 1 1 _1 j associated with the vortex lattice melting is 
clearly observed. Jl(J 

We first discuss the relaxation mechanism. It should 
be noted that the spin diffusion effect from the core re- 
gion is unimportant in the present case because of the 
large broadening of the spectra below Tm 10 1- Gener- 



ally, T]~ is composed of a contribution from vortex 
45 K. This result reinforces the fact that the observed vibration 17 T^y V and a contribution from the QP scat- 
frequency dependence of 17 T^ pl at lower temperatures tering 17 T^ 1 p 10 2lf|- 



originates from the field distribution associated with the 
vortex lattice. Figure 3 displays the T-dependencies of 
For the planar sites, we show the 
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largest and the smallest values of 17 T 1 pl for correspond- 
ing temperatures. The T-dependence of the 



The vortex vibration influences T\ through the trans- 
verse fluctuation of the local magnetic field. It has been 
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suggested that the vortex vibration plays an important 
role in T\ for very anisotropic Tl 2 Ba2Cu06 0]. We 
consider that in YBa 2 Cu40s, the QP scattering predom- 
inates 17 Ti at the planar sites, while the vortex vibration 
seemingly influences 17 Ti at the apical sites, based on the 
following explanation. First, as shown in Figs. 1 (a) and 
(b), l7 T{~ ap shows much smaller frequency dependence 
than 17 T^ pV Second, if the contribution from the vortex 
vibration dominates, then 17 should be the same as 
17 T 1 ~p i . However, as shown in Fig. 3, 17 T^l p are smaller 
than 17 T^ pl by a factor of ~ 5. Third, the vortex lattice 
melting seriously influences 17 Tj~ , while its influence on 
17 T- pl is fairly small. @ 

Having established that the 17 T^~ pl is dominated by the 
QP excitation, the next question is the frequency depen- 
dence. We first discuss l7 T^~ in the low frequency region 
well outside the core. Within the framework of the semi- 
classical approximation, it has been pointed out that the 
local DOS N(E — 0, r) well outside the core is domi- 
nantly produced by the Doppler shift of the delocalized 
QP spectrum, E(p) — > E(p) — v s ■ p, in the presence of 
the supercurrent flow with velocity | v s \= h/2mr. As 
a result, there is a finite N(E = 0,r)(oc| v s |), which 
decays as 1 jr for a single vortex and varies slightly more 
modestly than 1/r in the vortex lattice state. According 
to Ref. Q), this local DOS gives rise to a monotonous 
increase of 17 T^ pl with frequency. 

In the intermediate frequency regime where 17 T^~ pl in- 
creases, the Doppler shift appears to be important for 
the local DOS. However, the observed decrease of 17 T^~ pl 
when going from C- to B-point is apparently inconsis- 
tent with the Doppler shift. There are two possible ori- 
gins for this unusual behavior in 17 T^~ pl . One is the 
competition between the thermally excited QPs and the 
Dopplcr-shifted QPs in the presence of the strong an- 
tiferromagnetic fluctuation which is strongly suppressed 
by the Doppler shift, as pointed out in Ref. The 
other is the vortex lattice effect, which is not taken into 
account in the semi-classical theory of Doppler shift. Ac- 
cording to the numerical calculation, the localized QPs 
in the Andreev bound state around the vortex core ex- 
tend to (±7r, ±7r) directions, and can overlap the QPs of 
the neighboring vortices J22| . This effect can increase the 
QP DOS at the C-point, depending on the vortex lattice 
structure. At present it is not clear which effect is rele- 
vant to 17 T^~ pl at low frequency. The fact that the low 

frequency behavior of 17 T^ pl cannot be simply explained 
by the Doppler shifted QPs implies some key features of 
the QP excitation outside the vortex core. 

Finally we discuss 17 T^~ pl at high frequencies, at which 
17 T^~ pl begins to decrease after reaching maximum. It is 
apparent that this behavior is again incompatible with 
the Doppler shift. As pointed out in Ref. |l4j], the semi- 
classical description of the QPs breaks down in the vicin- 



ity of the core region (r < 2£) where the superconducting 
order parameter varies rapidly. Therefore it is natural to 
consider that the reduction of 17 T^ pl occurs in the vor- 
tex core regime. It should be noted that the antiferro- 
magnetic fluctuation, if it exists, always increases 17 T^ pl . 
These considerations lead us conclude that the decrease 
of 17 T^ pl is caused by the suppression of the QP DOS 
within the vortex core. This is the first experimental evi- 
dence of the reduction of the local DOS within the vortex 
core, which was probed by a bulk measurement. In the 
pure d-wave superconductors, the local DOS at the center 
of the vortex core should be enhanced, as it is in s-wave 
superconductors, owing to the interference of the Andreef 
scattered QPs at the both sides of the pair potential. It 
would follow that 17 T^ pl should be enhanced in the vor- 
tex core regime, which is apparently inconsistent with the 
present results. The present results are seemingly con- 
sistent with the STM measurements, which show that 
only one QP energy level is present within the core for 
the optimally doped YBa 2 Cu30y and Bi 2 Sr 2 CaCu 2 08+,s 
[EJ. Several possible origins for the reduction of DOS 
have been proposed as follows: (1) different symmetry, 
such as d x 2_ y 2 + id xy , is induced within the core ||; (2) 
the Andreef bound state is split due to the large A 2 jep 
value where A is the superconducting gap and ef is the 
Fermi energy; (3) the strong antiferromagnetic interac- 
tion changes the core state dramatically, as discussed by 
the calculation based on the t — J model [Q and the 
SO (5) theory [|. To clarify the origin of the reduction 
of the local DOS within the vortex core, more detailed 
information of the QP structure, especially in the under- 
doped regime, is strongly required. 

In summary, we measured the site-selective spin-lattice 
NMR relaxation rate in and around the vortex core of un- 
derdoped YBa 2 Cu408- Nonmonotonic NMR frequency 
dependence of T^ pl indicates that T^ pl well outside the 
vortex core regime cannot be simply explained by the 
Dopplcr-shifted delocalized quasiparticles of the d-wave 
superconductor and that the local DOS in the vortex core 
regime is strongly reduced. 
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Note added in proof. After submission (at April 12, 
2001) of this manuscript, we are aware of an article (in 
Nature 413 (October 4, 2001) 501) of similar NMR study 
on the spatial dependence of T^ 1 in ncar-optimally doped 
YBa 2 Cu 3 7 . 
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